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Abstract. In SUSY scenarios with invisible LSP, sparticle masses can be determined from fits to 
the endpoints of invariant mass distributions. Here we discuss possible improvements by using the 
shapes of the distributions. Positive results are found for multiple-minima situations and for mass 
regions where the endpoints do not contain sufficient information to obtain the masses. 
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INTRODUCTION 

In R-parity conserving supersymmetric models sparticles are produced in pairs at the 
collision point, then decay in cascades, resulting in a number of Standard Model particles 
as well as two Lightest Supersymmetric Particles (LSPs), one for each primary sparticle. 
If the LSPs leave the detector without a trace, as is the case in most scenarios of this 
type, the event cannot be fully reconstructed. This in turn prevents a direct measurement 
of the sparticle masses from mass peaks. The foreseen way to obtain these masses 
in such scenarios, is through endpoint measurements of invariant mass distributions 
. Particularly suited is the cascade decay qi — > — > Ijtlnq — » Xikkq 
which for a large share of SUSY models and parameters is kinematically allowed and 
appears at a rate sufficient for study. Four invariant mass distributions can be constructed 
from the visible decay products, m//, m q u, m q ^ high ^ and m q in QV/ ^. 

Once the endpoints have been measured, the masses can be found. However, measur- 
ing an endpoint is not trivial. A fit must be made to the edge region of the distribution. 
Such a fit necessarily involves a signal and a background hypothesis. If these are not in 
good correspondence with the true shapes, it is expected that the endpoint fitting may in- 
troduce large systematic uncertainties. While a straight-line hypothesis sometimesper- 
forms well, the fact that each distribution can take on a large variety of shapes o, [ZD 
really is a limiting factor, especially since a mismeasurement of an endpoint by 10 GeV 
can easily multiply to 40 GeV in the resulting masses. Without an appropriate signal hy- 



1 Additional constrained mass distributions can be constructed, e.g. m q u for events with m« > m™ ax / \/2. 
The corresponding endpoint measurements however usually involve large errors. Although relevant, such 
distributions have not been considered in this study. 



pothesis the endpoint method is incomplete. The many studies undertaken so far reflect 
this in that they have mainly focused on determining the statistical error, which can be 
fairly well estimated by a straight line, putting off the systematics of the fitting procedure 
for later investigations. With the recent advent of an analytic description of the shapes 
of these mass distributions BSD, the time for such investigations has come closer. 

MASSES FROM SHAPES 

The completion of the endpoint method requires formulas for the edge-region part of 
the distributions expressed by the endpoints. The shape formulas are however written 
in terms of the sparticle masses. Fitting the distributions gives the masses directly 
without any recourse to the endpoints. The considerable difference in approach asks for a 
comparison between the two methods. We here investigate the matter for the mSUGRA 
point SPS la [9]. In order for the characteristics of the methods to easily shine through, 
we make a simple comparison, with no background included, and no detector effects. 
The comparison is based on histograms generated randomly from the shape formulas 
taken from [SD convoluted with a gaussian of width 2 GeV for m\i and 10 GeV for the 
other distributions. The selected procedure and numbers mimic the effect of sparticle 
widths, initial and final-state radiation, detector effects etc. The same procedure is also 
invoked for the shape-fitting process. The normalisation of the histograms corresponds 
roughly to the ultimate statistics expected for an LHC experiment. 



TABLE 1. Endpoint values and errors [in GeV] for 
SPS la. 





Nominal 


Fitted 


Error 


Inflated error 


max 
m ll 


77.07 


77.13 


0.04 


0.04 


m max 


375.8 


378.9 


0.4 


1.2 




298.5 


304.2 


0.8 


2.4 


m max 
m qll 


425.9 


432.3 


0.3 


0.9 



Endpoint method. For the endpoint analysis a straight-line fit is used to obtain the 
statistical uncertainty. Table \T\ shows the nominal values together with the fitted values 
and the errors obtained. A second set of errors, three times the statistical ones, except 
for m™ ax , is included to give a slightly more realistic situation. The discrepancy between 
the nominal and the fitted endpoints reminds us of the problem with endpoint fitting. 

Figure \T\ (left and middle) shows the resulting 1, 2 and 3 a contours in the m^o-mj R 
plane for the two sets of errors. The nominal endpoint values are used, not the fitted 
ones. From these figures the main characteristics of the endpoint method are apparent: 
the occurrence of multiple solutions o, LL0D and strong mass correlation reflecting the 
fact that mass differences are accurately determined while the overall scale only poorly 
so. 

Shape method. Since the background is expected to be more prominent for lower 
invariant masses, only the higher part of the distributions are used in the shape fit. 
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FIGURE 1. SPS la comparison. Left (Middle): endpoint method for the first (second) set or errors. Two 
minima are found, the correct one situated at m^o = 96 GeV, the false at m^o = 80 GeV. Right: shape fit 
to the same distributions. The false minimum is not present. 



The result of the fit is given by the ^-function in Fig. Q] (right). Notice that the false 
solution is absent. The shape difference for the two mass sets is apparently sufficient 
to discriminate between them. Whether or not this situation holds for a more realistic 
analysis remains to be seen. The shape difference between these two points is below 
mSUGRA average The precision on the overall mass scale returned by the shape 
fit is roughly as for the endpoint approach using the second error set, except for the false 
solution present in the endpoint case which stretches to very low masses. Finally it is 
found that also for the shape fitting method the resulting masses are very correlated, 
constraining mass differences much more than masses. More study is however needed 
to compare the degree of correlation to that obtained from the endpoint method. 



WHEN ENDPOINTS ARE NOT ENOUGH 

In some regions of mass space the four endpoints are no longer linearly independent 
due to (w^ x ) 2 = (m™ ax ) 2 + ( m ™h igh )) 2 - When this happens, which is for regions (2,3), 
(3,1) and (3,2) in the notation of yD, curves exist in the four-dimensional mass space on 
which all mass points produce the same endpoints. In region (3,2) the curve is given by 

2 Xi 2 2 . \ m ll ) 2 2 . [m qll > A / V high 

m~ = — , mtn = m~ H , m~ = m~ H - , A = — 

i R 2 -A' *° A-l ' qL '* A-l ' /n™ 

(1) 

From one mass set new mass sets on the curve can be generated by keeping the end- 
points fixed, then changing m^o and calculating the other masses. Under these circum- 
stances four endpoint measurements do not impose sufficient constraints to obtain the 
masses. Extra information is however available from the shapes. Figure [2] (three left- 
most) shows the relevant gaussian-convoluted (width 10 GeV) mass distributions for the 
mass set (m^o,m^,m^o,m^ L ) = (100,300,500,600) GeV (black, solid), and two sets 
on the line defined by increasing (blue, dotted) and decreasing (red, dashed) m ? o 
by 30 GeV. While the endpoints are the same, the shapes are seen to differ. (The mu 
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FIGURE 2. Left to right: m q in^^\, ?m 9 /(i ow ) an d m q u for three points in region (3,2), all having the 
same endpoints, but different distributions. Right: % 2 -function from fitting the shapes has minimum at the 
correct m=o. Only smearing is applied, no statistical fluctuations, therefore Xmm = 0- 



distribution is not shown as its shape is always the same.) 

Based on the distributions of the original mass set (m^o = 100 GeV, black solid), his- 
tograms are generated, assuming a similar number of events as for the previous SPS la 
investigations. Figure [2] (right) shows the % 2 function from a shape fit to these distri- 
butions. The histograms contain no statistical fluctuations, which is why the minimum 
is at % 2 = 0. Adding statistical fluctuations will lift %^ from zero, but the conclusion 
stays the same: A clear minimum is found for the correct masses in this region where 
the endpoints alone can not bring the masses. 



CONCLUSION 

We have investigated the use of shape formulas to obtain the sparticle masses in sce- 
narios where the LSP is undetected. Comparison with the standard endpoint method 
was made for the mSUGRA point SPS la. The extra solution at m 9 o = 80 GeV which 
is returned by the endpoint inversion, is absent in the shape- fitting results. The strong 
correlation in the resulting masses remains: Mass differences are accurately determined 
while the overall mass scale has larger uncertainties. It was furthermore shown that the 
use of shapes allows for masses to be determined in regions of mass space where the 
endpoints alone do not suffice due to linear interdependence. In conclusion, the use of 
shapes to determine masses constitute a very promising approach. More study is needed. 
It is particularly important to understand what impact the distortion of the invariant mass 
distributions, particularly from selection cuts and various backgrounds, will have on the 
proposed shape fitting technique. 
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